THE BARNES G FUNCTION AND ITS RELATIONS WITH 
SUMS AND PRODUCTS OF GENERALIZED GAMMA 
CONVOLUTION VARIABLES 



A. NIKEGHBALI AND MARC YOR 

Abstract. We give a probabilistic interpretation for the Barnes G- 
function which appears in random matrix theory and in analytic num- 
ber theory in the important moments conjecture due to Keating-Snaith 
for the Riemann zeta function, via the analogy with the characteristic 
polynomial of random unitary matrices. We show that the Mellin trans- 
form of the characteristic polynomial of random unitary matrices and 
the Barnes G-function are intimately related with products and sums of 
gamma, beta and log-gamma variables. In particular, we show that the 
law of the modulus of the characteristic polynomial of random unitary 
matrices can be expressed with the help of products of gamma or beta 
variables, and that the reciprocal of the Barnes G-function has a Levy- 
Khintchin type representation. These results lead us to introduce the 
so called generalized gamma convolution variables. 



1. Introduction, motivation and main results 

The Barnes G-function, which was first introduced by Barnes in [3] (see 
also [I]), may be defined via its infinite product representation: 



G(l + z) = (2vr) 2/2 exp 



cxp 



z 2 



\[{i+i)z*+z\ n(i+ 

J n=l 

(1.1) 

where 7 is the Euler constant. 

From (II. ID . one can easily deduce the following (useful) development of 
the logarithm of G (1 + z) for \z\ < 1: 

log G (1 + z) = I (log (2tt) -l)_-(l + 7 )*a + £ (-I)™" 1 C {n - 1) £- 

n=3 

(1.2) 

where £ denotes the Riemann zeta function. 

This Barnes G-function has recently occurred in the work of Keating and 
Snaith [Hj in their celebrated moments conjecture for the Riemann zeta 
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function. More precisely, they consider the set of unitary matrices of size 
N, endowed with the Haar probability measure, and they prove the following 
results: 

Proposition 1.1 (Keating-Snaith [14]). If Z denotes the characteristic poly- 
nomial of a generic random unitary matrix, considered at any point of the 
unit circle (for example 1 ), then the following hold: 

(1) For A any complex number satisfying 9te(A) > —1: 



E 



N 



\z 



2A 



N 

n 



r (j) r (j + 2A) 
(r(j + A)) 2 



(1.3) 



(2) For^t(X) > -1; 



lim 



N x2 



E 



N 



\z 



2A 



(g(l + A)) 2 
G(1 + 2A) 



(1.4) 



Then, using a random matrix analogy (now called the "Keating-Snaith 
philosophy"), they make the following conjecture for the moments of the 
Riemann zeta function (see |14j.|16j): 



lim -5- — 

T^oo (l 0g T) A T 



dt 



2A 



M (A) A (A) 



where M (A) is the "random matrix factor" 

(G(l + A)) 2 



M (A) 

and A (A) is the arithmetic factor 
A (A) 



G(1 + 2A) 



n 



1 

1 - - 

p 




r(A + m) 
m!r(A) 



where, as usual, "P is the set of prime numbers. 

Due to the importance of this conjecture, as discussed in several papers 
in [16] , it seems interesting to obtain probabilistic interpretations of the non 
arithmetic part of the conjecture. More precisely, the aim of this paper is 
twofold: 

• to give a probabilistic interpretation of the " random matrix factor" 
M (A), and more generally of the Barnes G- function; 

• to understand better the nature of the limit theorem (jl.4p and its 
relations with (generalized) gamma variables. 

To this end, we first give a probabilistic translation in Theorem 11.21 of the 
infinite product (II. ip in terms of a limiting distribution involving gamma 
variables (we note that, concerning the Gamma function, similar translations 
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have been presented in [10] and [9] ) . Let us recall that a gamma variable 7(l 
with parameter a > is distributed as: 

F( 7a£ d t ) = t '"' e ^ ( r t)dt - «>° < L5 > 

r(o) 

and has Laplace transform 

E[exp(-A 7a )]= T^ra , «e(A)>-l (1.6) 
and Mellin transform: 

E[(7a)S] = £ r^T' *«(*)> -° ( L7 ) 

Theorem 1.2. // (7n) n >i ore independent gamma random variables with 
respective parameters n, then for z such that !SRe(z) > — 1: 

V 2N N "I 



lim — ^rE 



(1.8) 

where 

A = J±. (1.9) 

The next theorem gives an identity in law for the characteristic polynomial 
which shall lead to a probabilistic interpretation of the "random matrix 
factor" : 

Theorem 1.3. Let A denote the generic matrix ofU (N), the set of unitary 
matrices, fitted with the Haar probability measure, and Z]y (A) = det (I — A). 
Then the following hold: 

(1) For SHc(t) > —1, we have: 

j=i l± L7 ~r 2)) 

(2) Equivalently, in probabilistic terms: 

N N 

n^- =" i^(A)inv^ ( Ln ) 

i=i i=i 

where all variables in sight are assumed to be independent, and jj, 7'- are 
gamma random variables with parameter j. 

The Barnes G-function now comes into the picture via the following limit 
results: 

Theorem 1.4. Let ( r J n ) n >i be independent gamma random variables with 
respective parameters n; then the following hold: 



-1 
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(1) for any A, with £He(A) > — 1, we have: 



lim 




cxp 



-a$>(j)) =(A A G(1 + A) 



-1 



(1.12) 



(2) consequently, from ji.ig]) . together with hi. 11}) , we recover the limit 
theorem |i.^[ ) of Keating and Snaith: for fKe(A) > —1: 

n2 



1 



lim 



E 



TV 



2A 



(g(i+A)r 

G(1 + 2A) 



We then naturally extend Theorem 11,21 to the more general case of sums 
of the form: 



S N = Y,(Yn-^[Yn}) 



(1.13) 



71=1 



where 



1 / (n) 

2/i + 



Vn I 



and where (y, 



are independent, with the same distribution as 

l<j<n<oo 

a given random variable y, where Y is a generalized gamma convolution 
variable (in short GGC), that is an infinitely divisible K + -valued random 
variable whose Levy measure is of the form 



v (dx) 



dx\ 
x ) 



\i (d£) exp 



(1.14) 



where [i (d£) is a Radon measure on R + , called the Thorin measure associ- 
ated to Y . We shall further assume that 



which, as we shall see is equivalent to the existence of a second moment for 
Y. 

The GGC variables have been studied by Thorin [19] and Bondesson [6], 
see, e.g., [12] for a recent survey of this topic. 

Theorem 1.5. Let Y be a GGC variable, and let (Sn) as in tl.l3\) . We 

note 

a 2 = E[y 2 ] - {E[Y}) 2 . 
Then the following limit theorem for (Sat) holds: if X > 0, 

(1-15) 



N— 



■E [exp {-\S N )\ 



N^oo 



W(A), 
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where the function Ti (A) is given by: 
X 2 a 2 



H (A) = exp |^ (1 + 7 ) + (y) (exp (-Ay) - 1 + Ay - } 



with 



S M (y)= / M (d0 



2 sinh 



2 ' 



2 

(1.16) 
(1.17) 



Limit results such as fjl . 15j) are not standard in Probability theory, and 
we intend to develop a systematic study in a forthcoming paper [11]. The 
rest of the paper is organized as follows: 

• in Section 2, we prove Theorems 11.21 11-31 and 11.41 We also give an 
interpretation of Theorem 11.21 in terms of Bessel processes as well as 
an important Levy-Khintchine type representation for 1/G (1 + z). 

• in Section 3, we compare Theorems 11.21 and 11.41 which help us find 
adequate extensions for each of them (in particular we shall introduce 
there generalized gamma convolution variables and prove Theorem 
O); 



2. Proofs of Theorems 11.21 11-31 and 11.41 and additional 
probabilistic aspects of the barnes g-function 

2.1. Proof of Theorem 11.21 and interpretation in terms of Bessel 
processes. 

2.1.1. Proof of Theorem \l.S\ To prove Theorem ll.21 we simply use the fact 
that 



E 



exp 



,2™ 

n 



1 



1 + 



z\ n ' 



n 



Hence, for z > 0, the quantity 



E 



v 



exp 



\n=i 



equals 



We then write 
Dn exp 



n 



H — ) exp (— z) 




z 



v 2 N 



log N ) = Dm exp 



^ n 



exp 



logiV 



?1=1 




n=l 



which from (|l.lj) converges, as N — > oo, towards 



G(l + z) (2^)" 2/2 exp 



z + z' 
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which proves formula (|1.8p . for z > 0, and thus for *He(z) > — 1 by analytic 
continuation. 

2.1.2. An interpretation of Theorem in terms of Bessel processes. Let 
(i?2n (*))*>o denote a BES(2n) process, starting from 0, with dimension 2n; 
we need to consider the sequence (i?2n) n=1 2 °f such independent processes. 
It is well known (see [17] for example) that: 



for fixed t > 0, R\ n (t) 



law 



(2*)7n- 



Moreover, we have: 



for fixed t > 0, R 2 2n (t) = 2 j ^/i? 2 n (s)d^ n) + 2nt 

the stochastic differential equation of -R 2 „, driven by a Brownian Motion 
Thus 



N 



E R 2n CO If i [ 7n 



n=l \n=l 

We now write Theorem 11.21 for 9\t(z) > —1/t as: 



lim 

N^oo 



N— 



-E 



exp — z 



E 

71=1 



Rln (0 " 2nt 



2n 



A exp 



where 



G(l + tz) 
(2.1) 



V 2?r 

We now wish to write the LHS of (|2.ip in terms of functional of a sum of 
squared Orstein-Uhlenbeck processes; indeed, if we write: 



N 

E 

n=l 



Rln (*) ~ 2nt 
2ra 



2 Ei ^^Lw^»>, 



the RHS appears as a martingale in i, with increasing process 



N rt 1 



(s)ds, 



and we obtain that (|2.1|) may be written as: 



lim 



1 



6XP ( _ |S / i^in( s ) ds 



^4 exp 
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where, under P* 2 ) the process (i?2n (*))t>o satisfies (by Girsanov's theorem) 




= 2 f Ji?| n ( S )d^) - - / R\ n (s) ds + 2nt 
Jo n Jo 

That is, under p(*), (i?| n (i)) t>0 now appears as the square of a one dimen- 



sional Ornstein-Uhlenbeck process, with parameter 




2.2. Proof of Theorem [L3l Formula (fTTTUj) is due to Keating-Snaith [H] . 
Formula (jl.lip follows from (jl.lOp once one recalls formula (jl.7p . 

2.3. The characteristic polynomial and beta variables. One can use 

the beta-gamma algebra (see, e.g., Chaumont-Yor [8], p. 93-98, and the refer- 
ences therein) and (jl.lip to represent (in law) the characteristic polynomial 
as products of beta variables. More precisely, 

Theorem 2.1. With the notations of Theorem \1.3[ we have: 

i^i 's- 2 « (nv^i) (n v^iw) < 2 - 2 > 

where all the variables in sight are assumed to be independent beta variables. 

Proof. To deduce the result stated in the theorem from (jl.lip . we use the 
following factorization: 

7, = W yf^fr, (2-3) 

where on the RHS, we have: 



for j = 1 
for j > 1 




Indeed, starting from (|2,3p . then multiplying both sides by yljl'j and using 
the beta-gamma algebra, we obtain that (|2.3|) is equivalent to: 

law / ; — 

which easily follows from the duplication formula for the gamma function 
(again, see, e.g., Chaumont-Yor [8], p. 93-98). □ 

Remark 2.2. In the above Theorem, the factor 2^ can be explained by the 
fact that the characteristic polynomial of a unitary matrix is, in modulus, 
smaller than 2^. Hence the products of beta variables appearing on the 
RHS of the formula (j2.2p measures how the modulus of the characteristic 
polynomial deviates from the largest value it can take. 
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2.4. A Levy-Khintchine type representation of \/G{l + z). In this 
subsection, we give a Levy-Khintchine representation type formula for 1/G (1 + z) 
which will be used next to prove Theorem II. 41 



Proposition 2.3. For any z £ C, such that $tt(z) > —1, one has: 

i J2 poo du I exp (—zu) — 1 + zu 



i z z r 

exp ) __ (i g (2vr) - 1) z + (1 + 7 ) — + jf 



+ ^1 2 ^ toV "' v ^~ 1 u 2 j Q u(2sinh(u/2)) 2 

(2.4) 

Before proving Proposition 12.31 a few remarks are in order. 



Remark 2.4. Note that formula (|2.4p cannot be considered exactly as a 
Levy-Khintchine representation. Indeed, the integral featured in (|2.4p con- 
sists in integrating the function ^exp (—zu) — 1 + zu — ^-^J against the 

measure ^ , which is not a Levy measure. Indeed, Levy mea- 

u(2sinh(tt/2)f 

sures integrate (u 2 Al), which is not the case here because of the equivalence: 
u (2 sinh (u/2)) 2 ~ u 3 , when u — > 0. Also, due to this singularity, one cannot 
integrate (exp (—zu) — 1 + zul^ u< ^) with respect to this measure, and one 

2 2 

is forced to "bring in" under the integral sign the companion term '^ L -^. 
Proof. From the consideration of the series 

oo 



C(z):=^(-l) n ^((n-l)- (2.5) 

n=3 



featured in (jl.2p . it seems natural to introduce a random variable Q taking 
values in M + , with Mellin transform: 



ER'lHV r^fM (2-6) 




where e denotes a standard exponential variable. A little more analytically, 
formula ()2.6[) may be presented as 

E[Q S ] = Ar( s + 3)C(s + 2). (2.7) 

We first show the existence of the random variable Q by computing its 
density (for an example of occurrence of the random variable Q in the theory 
of stochastic processes and some relation with the theory of the Riemann 
Zeta function, see e.g. [5] and [1]). 

From the definition of Q via its Mellin transform, we get for / : M + — > M_|_: 
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7T 2 



vn=l 



Hence by some elementary change of variables: 

n exp (—nu) du. 

1 

Now, since Y2ri=l ex P (~ nu ) 

Y.n=i nex V (~ nu ) = ~, 7 

(2smh lf/ )j 

u/2 



1 

exp (u) — 1 



, we get: 



Hence, 



1 (Q G du) 



7T 



sinh(|) 



du 



(2.8 



Now, we consider the following series development for \z\ < 1: 



E 



Q 3 I ^ 



vn=l 



/?! 



4EH) n c(«-i)-. 



n=3 



hence, in comparison with formula (|2.5p we obtain: 



E 

3 



Q 3 



n=l 



71! 



vr 2 

E 

3 



—3 ( exp {-zQ) -1 + zQ 



Q 2 z 2 



Now, formula (|2.4p . for |z| < 1, follows from (12. 8p and the fact that: 
1 



G(l + z) 



exp(logG(l + z)) 

z 



exp 



1 



2 (log(27r)-l) + ^(l+ 7 ) z 2 -C(z) 



The formula extends by analytic continuation to the case 9^e(z) > — 1. □ 

2.5. Proof of Theorem 11.41 To prove Theorem 11.41 we shall use the fol- 
lowing lemma: 

Lemma 2.5. For any a > 0, the random variable log (7 tt ) is infinitely di- 
visible and its Levy-Khintchine representation is given, for SHe(A) > —a, 
by: 

r (a + A) 



E 



where 



T(a) 



exp < \if) 



r 

■(«)+ / 



exp (—an) (exp (—An) — 1 + An) 



n (1 — exp (—it)) 



du 



r(a)- 
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Proof. This is classical: it follows from (11.70 and some integral represen- 
tation of the ^-function, see, e.g. Lebedev [15] and Carmona-Petit-Yor [7] 
where this lemma is also used. □ 

We are now in a position to prove Theorem 11.41 We start by proving the 
first part, i.e. formula (|1.12p . Let us write 

A" 

/ -v \ 

MA):= ' 



N x2 / 2 



E 



N 



Thus with the help of Lemma 12.51 we obtain: 



A V V (J) + / ( { =1 ' (exp (-An) - 1 + An) du 

j^y Jo u(l-exp(-u)) 



Note that: 



(2.9) 



N 



ll 



(1 - exp (-«)) 



exp (— u) (1 — exp (—Nu)) 
u (1 — exp (— u)) 2 



and we may now write: 



N 



where 
Jat(A) 



In (A) = exp < 



duexp (— u) 



3=1 



(2.10) 



2 0- - exp (-Nu)) (exp (-Xu) - 1 + Xu) } -— log N 

o u(l — exp(— u)) 



A2 
2 



Next, we shall show that: 



(2.11) 



JnW 



Joo(X) 



together with some integral expression for (A), from which it will be 
easily deduced how (A) and G (1 + A) are related thanks to Proposition 



We now write (|2.1ip in the form: 
duexp (— u) 



Jn(X) 



X 2 u 2 

2 v 1 - ex P (~Nu)) ( exp (-Xu) -1 + Xu 

o u (1 — exp (—u)) V l 



X 2 



duu exp (—u) 



2(l-exp(-iVn))-logiV 

2 \Jo (1- exp (-it)) 



Now letting N — » oo, we obtain: 

■Tiv (A) — > Joo (A) , 



(2.12) 
(2.13) 
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with 

Joo(A) = / — y \ —2 I exp (—Xu) — 1 + Xu — J + —C 

Jo u(l — exp(—u)) \ 2)2 

= Joo(A) + yC 



where 



N->oo \J (i _ exp (-it)) 



a limit which we shall show to exist and identify to be 1 + j(= C) with the 
following lemma: 

Lemma 2.6. We have: 

WeXp( ~ K) (1 - exp (-Nu)) du = log N + (1 + 7 ) + o(l). 
Jo (1 — exp (—it)) 

Consequently, we have: 

C= l+ 7 . 

Proof. 

JV „oo AT „oo oo 

= > / - r- exp (— ku) du = > / cfaz exp ( — fett) u > exp(— ru) 

(l-exp(-u)) ^Jo h 

N 

= 22/ ^uexp(-(r + fc)it)=^^ 2 

7^ nJo iT n\r + k) 

k=l r=0 u fc=l r=0 v ' 

Now, we write: 



N oo 1 AT-1 / k "> 

AT-1 fc 



AT-1 fc 

.s- 



fc=l s=l 



S 2 
s=l 

Ar A 1 1 " 1 1 

^ s 2 N ^ s N 

s=N+l s=l 
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The result now follows easily from the facts: 

1 



hm JV V ^ 

Af^oo S 2 
s=N+l 

N 



logiV + 7 + o(l) 



□ 



We have thus proved so far that: 



lim 



1 



N^oo N XZ I 2 

where 

Joo(A) 



E 



N 



N 



exp [ -Aj^V (j) ] = exp ( (A)) , (2.14) 



ditexp (— u) 
/o it(l — exp (— u)) 
We can still rewrite Jqo (A) as: 

Joo(A) ' 



exp (—Ait) — 1 + Ait 



AV\ a 



+ y (1 + 7)- 



^2 M 2\ ^2 

2 ( exp (-Ait) - 1 + Ait — ) + — (1 + 7) . 



la it(2sinh(ii/2)) 
Now comparing (|2.4p and (|2,15p . we obtain: 

exp(J oc (A)) = (A A G(1 + A) 
Plugging (|2.16p in (|2.14j) yields the first part of Theorem 11.4 



(2.15) 
(2.16) 



1 



lim , , , 



E 



N 



V 



exp [ -A ^ ^ (i) 



1 



e G(l + A) 



To prove the second part of Theorem 11.41 we use formula (jl.lip together 
with formula (|1.12j) . Formula (jl.lip yields: 

2n 



1 



E 



N 



2A 



1 



ATA 2 



E 



v 



J N\ 



|2A 



/ 



1 



E 



V" 



N 



(2-17) 

Multiplying both sides by exp ( — 2A X^fci V* 0) ) an d using (|1.12p we obtain: 



1 



lim 



E 



A' 



\z 



2A 



(G{l + \)f 
G(1 + 2A) 



which completes the proof of Theorem 11.41 
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3. Generalized Gamma Convolutions 
3.1. Definition and examples. We recall the definition of a GGC variable. 

Definition 3.1. A random variable Y is called a GGC variable if it is 
infinitely divisible with Levy measure v of the form: 

(IT f 

u(dx) = — Md£)exp(-£x), (3.1) 
x J 

where [i (d£) is a Radon measure on R + , called the Thorin measure of Y . 

Remark 3.2. Y is a selfdecomposable random variable because its Levy mea- 

dx 

sure can be written as v (dx) = — h(x) with h a decreasing function (see, 
e.g. [IB], p.95). 

Remark 3.3. We shall require Y to have finite first and second moments; 
these moments can be easily computed with the help of the Thorin measure 

M(d£): 



a 2 



E[Y] = 1^-1 = J M(d£)~ 
E[Y 2 ]-(E[F]) 2 = /x_ 2 = |/x(dC)^ 



Now we give some examples of GGC variables. Of course, j a falls into 
this category with \i (d£) = a^i (d£) where 5i (d^) is the Dirac measure at 1. 
More generally, the next proposition gives a large set of such variables: 

Proposition 3.4. Let f be a nonnegative Borel function such that 

\ du log (1 + /(«)) <oo, 
Jo 

and let (7^) denote the standard gamma process. Then the variable Y defined 
as 



Y= / d 7tt /(«) (3.2) 
J 

is a GGC variable. 

Proof. It is easily shown, by approximating / by simple functions that 

/ f°° f 00 dx \ 
E[exp(— AY)] = exp ( — J du J — exp (— x) (1 — exp (— A/ (u) x)) \ 

= exp ("i f(i d " exp ("7ck)) (1 " exp( " Ay)) 

which yields the result. □ 
For much more details on GGC variables, see 1121 . 
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3.2. Proof of Theorem II. 51 We now prove Theorem 11.51 which is a nat- 
ural extension for Theorem 11.21 Recall that in (|1.13p . we have defined SV 
as: 

N 

S N = J2(Xn-®[Yn}) 

n=l 

where 



n 

and where \yf^ ) are independent, with the same distribution as 

a given GGC variable Y, which has a second moment. For any A > 0, we 
have: 

E [exp (— AiSjy)] = J] {$ ) ( 3 - 3 ) 

where 

£(A)=E[exp(-A(F-E[Y]))]. (3.4) 

Now, using the form (|3.1|) of the Levy-Khintchine representation for Y, we 
obtain: 

n (*(»)) = e ^{-E n J^ d ^( i -^ x - e ^(~ x ))) 

= exp(- J n(dZ)I N (Z,\)\ . (3.5) 

where: 

N f°° dx fX ( A \ \ 

In(£A) = Xj n _/ — exp(-^x)^l - -x - exp xj J 

= I — ex P(-<2/) (!--%- ex P (--%)) 

ti Jo y 

f°° dy ( N \ 
= I — y]rcexp(-n£y) (1 - Ay - exp (-Ay)) 

Jo v \ti J 

Some elementary calculations yield: 

E. exp (a) (1 — exp (— aN)) Nexp(-aN) 
nexp(— na) = 5 ; ^ — . 
n=1 (exp(a)-l) 2 (exp(a)-l) 

Consequently, taking a = £y in the formula for Ijy (£j A), we can write it as: 
Iat(£,A) = J A r(C,A)-i? JV (e,A) (3.6) 
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where 
Jn (£, A) 

and 

Rn (£, A) 
It is clear that 



dy 

o y 



1- Xy- exp (-Ay)) 



exp (£y) (1 - exp (-£yN)) 



dy 
o y 



(1 - Ay - exp (-Ay)) 



lim i?Ar(C,A) = 0. 

iV— »oo 



(exp (£y) - l) 2 

iVexp (-£yN) 
(exp (fy) - 1) 



(3.7) 
(3.8) 



Now let us study Jjv (£, A) when N — > oo. Let us "bring in" the additional 
term — more precisely, we rewrite Jn (£, A) as: 



■M£,A) 



°° dy 

o y 



1 x a. A V r x ^ 

1 - Ay + — exp (-Ay) 



exp (gy) (1 - exp (-&/AQ) 
(exp (£y) - l) 2 



/ d y y^ P (1 " 6XP (-^) 



(exp (£y) - l) z 



(3.9) 



Hence: 

J iV (e,A)=o(l)+ 



dy 



1 



where 



^(e) 



o y (2sinh(£y/2)) 

exp (£y) (1 - exp (-£yiV)) 



A 2 y 2 \ A 2 

, | 1 - Ay + — exp (-Ay) ) — —K N (£) 



2 

(3.10) 



dy y- 



1 



(exp (£y) - 1 
uexp (—it) 



e Jo (1- exp (-«))' 



(1 - exp (-Nu)) du. 



Moreover, from Lemma 12.61 we have: 

1 



K N (0 = £2 VogN + 1 + 7 + o(l)) . 



(3.11) 



Now using the above asymptotics, and integrating with respect to // (d£) 
(see (|3.5p ). we obtain: 



N 

III 



n=l 



exp |o(l) + ^/i_ 2 (log AT + 1 + 7 ) + jH (y) f 1 - Ay + ^- - exp (-Ay)) | 



(3.12) 
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Hence we have finally proved that if A > 0, 

— ^ E [exp (-ASW)] ► H(A), (3.13) 

where the function 7i (A) is given by: 

U (A) = exp | ^f- (1 + 7 ) + (y) (exp (-Ay) - 1 + Ay - ^) } , 

(3.14) 

with 




(3.15) 



which is Theorem 11.51 
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